INFLUENCE OF SEDIMENTATION ON CONVECTIVE INSTABILITIES IN 

COLLOIDAL SUSPENSIONS 
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We investigate theoretically the bifurcation scenario for colloidal suspensions subject to a vertical 
temperature gradient taking into account the effect of sedimentation. In contrast to molecular 
binary mixtures, here the thermal relaxation time is much shorter than that for concentration 
fluctuations. This allows for differently prepared ground states, where a concentration profile due 
to sedimentation and/or the Soret effect has been established or not. This gives rise to different 
linear instability behaviors, which are manifest in the temporal evolution into the final, generally 
stationary convective state. In a certain range above a rather high barometric number there is a 
coexistence between the quiescent state and the stationary convective one, allowing for a hysteretic 
scenario. 

Keywords: binary fluids; Soret effect; barometric number; coexisting stationary solutions; hysteresis. 



1. Introduction 

Thermal convection in binary mixtures has attracted 
much research activity already for along time [Flatten 
& Legros 1984, Cross & Hohenberg 1993, Liicke et al. 
1998]. In comparison to the pure fluid case, the dynam- 
ics and the bifurcation scenario are more complicated due 
to the extra degree of freedom associated with the con- 
centration field. Thereby solutal currents are not only 
driven by concentration gradients, they occur also in re- 
sponse to temperature inhomogeneities. This is denoted 
as the thermo-diffusive or Soret effect. Its influence on 
the convective buoyancy force is quantifled by the di- 
mensionless separation ratio "0. The sign of ip indicates 
whether temperature- and solutal-induced density gradi- 
ents are parallel or opposed to each other. At negative ^ 
the motionless conductive state experiences an oscillatory 
instability, saturating in a nonlinear state of traveling 
waves [Liicke et al. 1998]. On the other hand, at positive 
ip the convective instability remains stationary, but the 
critical Rayleigh number for the onset of convection is 
dramatically reduced as compared to the pure-fluid ref- 
erence value i?a° = 1708. This is a result of the joint 
action of thermal and solutal buoyancy forces. 

A typical property of binary mixture convection is the 
formation of concentration boundary layers [Winkler & 
Kolodner 1992]. This is a consequence of the fact that 
the concentration diffusivity D,. in mixtures is usually 
much smaller than the heat diffusivity n. For molec- 
ular binary mixtures the dimensionless Lewis number 
L = Dc/n adopts typical values between 0.1 and 0.01 
[Kolodner 1988]. If colloidal suspensions are under con- 
sideration, the time scale separation is even more dra- 
matic. In this context magnetic colloids, known as fer- 
roffuids, are a canonical example. These materials are 
dispersions of heavy solid ferromagnetic grains suspended 
in a carrier liquid [Rosensweig 1985]. With a typical di- 
ameter of 10 nm the particles are pretty large on molecu- 
lar length scales, resulting in an extremely small particle 
mobility. This feature is reflected by Lewis numbers as 
small as i = 10~^ [Blums et al. 1997]. The smallness of 



L leads to a situation where de-mixing effects take place 
on very large time scales. Thus, in those experiments, 
where thermodiffusion is irrelevant, ferrofluids can safely 
be treated as single-component fluid systems. 

However, ferrofluids and other colloidal suspensions 
are also known to exhibit a very large separation ratio 
tp (up to IVj] « 100 [Blums et al. 1999, Lenglet et al. 
2002] . This observation is due to the pronounced thermo- 
diffusivity of these materials in combination with the 
large specific weight difference of the two constituents. 
As a result, in these materials the solutal buoyancy forces 
are rather strong and a two-component treatment of con- 
vective instabilities is mandatory. By considering the 
classical Rayleigh Benard setup it is shown [Ryskin et 
al. 2003] that the convective behavior is significantly 
different from the case of molecular mixtures. Start- 
ing from the motionless configuration with an initially 
uniform concentration distribution, convective perturba- 
tions are found to grow even at Rayleigh numbers well 
below the threshold i?a" of pure-fluid convection. The ac- 
tual critical Rayleigh number Rac is drastically smaller, 
but experimentally inaccessible due to the extremely slow 
growth of convection patterns for Ra > i?ac, requiring 
very large observation times. On the other hand, op- 
erating the colloidal convection experiment at Rayleigh 
numbers Rac < Ra ~ Ra^, reveals considerable positive 
growth rates, which lead to a saturated nonlinear state 
almost as fast as pure- fluid convection does at Ra > Ra'^. 

In an external magnetic fleld the apparent imperfection 
of the bifurcation is even more pronounced in the case 
of ferrofluids. Magnetophoretic effects as well as mag- 
netic stresses have been taken into account in the static 
and dynamic parts of the equations leading to rather pro- 
nounced boundary layer profiles (with respect to the con- 
centration and magnetic potential) . This boundary layer 
couples effectively to the bulk behavior due to the mag- 
netic boundary condition [Ryskin & Fleiner 2004]. 

In the case of a negative separation ratio (negative 
Soret coefficient) the thermal and solutal density gradi- 
ents are opposed when heating from below. The linear 
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convcctivc oscillatory instability known from molecular 
binary mixtures (with tjj < —1) at Ra^ is also found for 
colloidal ones, but the nonlinear treatment shows that 
the linearly unstable oscillatory states are transients only 
and decay after some time, rendering the final convection- 
free state stable [Ryskin & Pleiner 2005] . Above a second 
threshold, somewhat higher than Ra°, a finite amplitude 
stationary instability is found, while small amplitude dis- 
turbances do not destroy the convection-free state. The 
traveling wave solution dominating in molecular binary 
mixtures is shifted to unrealistically high temperature 
gradients and is not possible anymore in colloidal sys- 
tems [Huke et al. 2000, Huke et al. 2007]. When heating 
from above molecular binary mixtures with a negative 
separation ratio tp < —1, a. linear stationary instability 
is found, which is basically driven by the solutal buoy- 
ancy and only slightly modified by thermal variations. 
In colloidal suspensions, however, the concentration and 
temperature dynamics show completely different behav- 
ior. Thus, this stationary instability is very different from 
that obtained by heating from below with a positive sep- 
aration ratio. In the former case small scale structures 
arise at very high Ra numbers, whose wavelength de- 
creases strongly with increasing Ra. 

In earth's (vertical) gravity field, the density contrast 
in colloidal suspensions results in a tendency to phase 
separate the two constituents. However, for truly col- 
loidal systems the particles are small enough that Brow- 
nian motion successfully prohibits a real phase separa- 
tion allowing for the binary mixture description. (In the 
different case of micrometer-sized particles, e.g. magne- 
torheological fluids, a two-fluid description [Onuki 1989, 
Milner 1989, 1993, Pleiner & Harden 2003] should be 
used.) Nevertheless, there is a slight accumulation of 
the heavier constituent towards the bottom, i.e. a sed- 
imentation induced concentration gradient, that can be 
expected to be relevant in systems with a large separa- 
tion ratio ^, while in molecular binary mixtures this ef- 
fect generally is negligible. Some experiments show that 
sedimentation strongly affects the qualitative behavior 
of thermal instabilities in ferrofluids (Bozhko & Putin 
[2003], Bozhko et al. [2006], Tynjala et al. [2006]). A 
first theoretical discussion has been provided by Shliomis 
& Smorodin [2005]. 

In ferrofluids instead of the temperature gradient an 

external magnetic field produces a destabilizing force (the 
Kelvin force), which is larger in areas with a higher con- 
centration of magnetic particles. Concentration fluctua- 
tions are therefore amplified and can lead to an insta- 
bility. Recently, Ryskin & Pleiner [2007] have shown 
theoretically that a gravity-stratified ferrofluid indeed be- 
comes convectively unstable in a sufficiently strong exter- 
nal vertical magnetic field. The amplitude of the veloc- 
ity field is a rather complicated fimction of time. Initially 
the amplitudes grow exponentially with the linear growth 
rate. For all realistic parameter values the final flow at 
long times is stationary. For intermediate times the am- 
plitude saturates at a high value and then decreases con- 



siderably to its asymptotic value. The transition between 
the intermediate high peak value and the very small, sat- 
urated one is due to the fact that the convective flow ef- 
fectively reduces the concentration gradient, thus reduc- 
ing the very basis for the instability. The final station- 
ary state is reached, when the process of building up the 
concentration gradient due to sedimentation is balanced 
by its destruction due to advection. Since the former 
process is very slow, only a very small velocity is neces- 
sary. In contrast to the velocities the amplitudes of the 
concentration variations are not small. The stationary 
concentration profile is essentially neither linear nor con- 
stant due to the nonlinear part of the Kelvin force. This 
is in marked contrast to the thermal convection problem 
in non-magnetic colloidal suspension. 

In this paper we consider a horizontally infinite slab 
of colloidal suspension (thickness h) subject to gravity 
and a vertical temperature gradient (no magnetic field). 
We restrict ourselves to the case of a positive separation 
ratio and heating from below. We provide a comprehen- 
sive discussion of the influence of sedimentation on the 
bifurcation scenario for that case. Sedimentation affects 
convective instabilities in two stages. First, on the level of 
possible ground states and their linear stability, and sec- 
ond on the long-time nonlinear asymptotic . Due to the 
slow particle diffusion one can prepare different ground 
states. Applying the temperature gradient to a suspen- 
sion right after its preparation, sedimentation has no time 
to develop (" homogeneous state" ) and the linear instabil- 
ity behavior is the same as in the case without sedimenta- 
tion [Ryskin et al. 2003]. However, waiting long enough 
until the sedimentation has taken place (" stratified case" ) 
and applying then the temperature gradient, leads to sit- 
uation similar to that of a negative tp material (and dis- 
regarding sedimentation): the destabilizing temperature 
gradient has to overcome the stabilizing sedimentation- 
induced concentration gradient, which leads to an en- 
hanced threshold and an oscillatory linear instability. At 
an even stronger temperature gradient the oscillation fre- 
quency goes to zero rendering the linear instability to be 
stationary. (This latter feature has no counterpart in the 
case of negative tp without sedimentation.) After having 
set up the mathematical framework in Sec. 2, the linear 
instabilities will be discussed in more detail in Sec. 3. 
Here we also have included the experimentally hardly re- 
alizable case that the concentration gradient, due to sed- 
imentation as well as due to the thermodiffusion effect, 
is fully developed ("fully developed state"). All these 
linear cases influence the nonlinear bifurcation scenario 
and are manifest as transients. In Sec. 4 we give an 
approximate analytical solution for the stationary non- 
linear problem. We predict a hysteretic behavior due to 
the effects of sedimentation in a parameter range, where 
the non-convective state coexists with finite amplitude 
stationary convection state. Numerically we describe the 
transients to the stationary state, find an oscillatory so- 
lution, and give analytical conditions, for which such a 
solution can exist. These main results are summarized in 
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Sec. 5. 

2. Basic Equations 

We consider a slab of colloidal liquid subject to a positive 
temperature gradient (in z direction) and gravity (in -z 
direction). The system of equations is the same as for 
molecular binary mixtures including incompressibility, as 
well as momentum, heat, and mass conservation reading 
in Boussinesq approximation [Boussinesq 1903, Flatten 
& Chavepeyer 1976, Brand et al. 1984] 



V • V 

dfV + V ■ Wv 

dtT + v- VT 
dtC + v- VC 



0, (1) 
-VW + PrV^v 

+Pr Ra [{T - T) - ^(C - C)] (2) 

V^T, (3) 

^(V^C + V^T). (4) 



The relevant variables are the flow velocity v, temper- 
ature T, and concentration (of the particles) C. The 
material is characterized by the transport coefficients k, 
Dc, Ds, and v for heat and particle diffusion, thermod- 
iffusion, and viscosity, respectively. As usual they are 
used to make the equations dimensionless scaling length 
with the layer thickness h, time with h? / k, tempera- 
ture with AT, the applied temperature difference, and 
concentration with {Dg/Dcj^T. The quantities T and 
C are reference values defined as the mean values for 
temperature and concentration. The pressure is a 
kind of Lagrange multiplier that serves to guarantee in- 
compressibility for all times. We are left with three di- 
mensionless numbers governing the bulk material prop- 
erties, the Prandtl number Pr = v/k, the Lewis number 
L = Dc/k , and the separation ratio tp = Dsf3c/ (DcPt), 
where /3t = —{1/ p)dp/dT and (ic = {1/ p)dp/dc are the 
thermal and solutal expansion coefficients. The dimen- 
sionless Rayleigh number Ra = Pxgh^ / [kv) , with g 
the earth's gravity constant, is the control parameter for 
the bifurcation behavior. 

Gravity is not only responsible for the driving force, 
but also for an inhomogeneous distribution of the par- 
ticles. In equilibrium, the balance between gravity and 
Brownian motion leads to a Boltzmann distribution for 
the concentration [Biben et al., 1993] with the sedimen- 
tation length hg — ksT/ (xTiTipg), where is the effec- 
tive buoyant mass of a particle, and xt is the osmotic 
compressibility. Since usually hs ^ /i, the exponential 
distribution reduces to a linear concentration profile due 
to sedimentation 



Ciz)^Cil-z/h,), 



(5) 



where C is the mean mass fraction of the colloidal par- 
ticles. To reach this true equilibrium state one has to 
wait for a rather long time allowing experiments to be 
started from either this state or from the homogeneous 
state C = C. The transition from the former to the lat- 
ter state is due to a mass flux of particles, which can be 



written in the simplest (and dimensional) form as [Blums 
2002] 



Dr 



(6) 



This sedimentation current should be added to the dif- 
fusive and thermo-diffusive concentration currents in the 
concentration dynamics Eq. ^ , but is generally neglected 
there by putting G = C resulting in ^ • j s = 0. At the 
rigid (and impermeable) boundaries, however, the sedi- 
mentation current cannot be neglected w.r.t. the other 
concentration currents, since there the total concentra- 
tion current has to vanish, which is guaranteed by the 
(dimensionless) boundary conditions 



(a,c + a,T)u^±i/2 
■yU=±i/2 

T\z=±l/2 



B 



-0 Ra 



0, 



(7) 
(8) 
(9) 



at the upper (z = 1/2) and the lower {z = —1/2) plate. 
The additional standard boundary conditions for flow 
and temperature reflect the no-slip condition and the ex- 
ternally applied temperature gradient across the mean 
temperature T. The boundary condition for the con- 
centration variable contains the barometric number B 
[Shhomis & Smorodin 2005] 



B 



(10) 



which can vary considerably, typically from 1 to 10^, 
due to the strong /i-dependence. This parameter mea- 
sures the importance of sedimentation relative to viscos- 
ity and thermal diffusion. A strong impact of sedimen- 
tation on the bifurcation scenario can be expected for 
B > ij} Ra^ . For ferrofluids, using a magnetic field Hq 
instead of the temperature gradient, a magnetic baro- 

metric number Bm = Xc^ h'^ I {P^^'^^D governs the 
sedimentation effects on convective instabilities [Ryskin 
& Fleiner 2007]. Here, Xc = dx/dC describes the con- 
centration dependence of the magnetic susceptibility x 
giving rise to magnetophoresis, and e is the effective mag- 
netic permeability of the material. 

3. Linear Instabilities 

Starting from the homogeneous state (no sedimentation 
gradient) the linear development of a convective insta- 
bility is unaffected by sedimentation effects, since there 
is no concentration gradient to be adverted. Thus, the 
true threshold, Rac is by a factor of ij} smaller than in 
the single-component case, i?a|^, but the time evolution 
is too slow for the instability to be observed. At higher 
Ra numbers the growth rates are sufficiently high, but 
the system is already in the instability regime mimicking 
an imperfect bifurcation to a stationary state. 
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In the stratified state (with a sedimentation-induced 
concentration gradient) the instabiUty threshold is higher 
than in the case before, since the concentration gradient 
opposes the temperature gradient. Using the analogy to 
the case of binary mixtures with a negative separation 
ratio one can expect a linear oscillatory instability above 
RoP^. However, in contrast to the thermo-diffusive con- 
centration gradient, the sedimentation-induced one is not 
proportional to the temperature gradient. Thus if the 
barometric number is too small, sedimentation cannot 
compete with the temperature gradient and a stationary 
instability has to expected. Indeed a linear stability anal- 
ysis along the lines of Ryskin & Pleiner [2005] gives the 
threshold condition 

ma Pr ^ + = ^BPr + 27TJ:'^Pr\ + 7X^ (11) 

with one additional sedimentation contribution (~ B). 
Here, the growth rate A = and = ilu for the stationary 
and oscillatory instability, respectively. The former oc- 
curs for small barometric numbers B < Bi = LijjRa^ at 
the threshold 



B 



(12) 



independent of any wave number, while for _B > _Bi an 
oscillatory linear instability is found with 



Ra° 



Ra" 



27Pr 



14 + 21Pr 



-B 



QPr 



14 -I- 27 Pr 



[B-B,] 



(13) 
(14) 



where the critical wave number approximately hy kc = 
and i?a" = IStt^. For typical material parameters, 
-0 ~ 10 and L ^ 10~^, the transition occurs at Bi ~ 1. 
Equations (fT3)) and p4)) were (for ■0 = 0) also obtained 
by Shliomis & Smorodin [2005], with slightly different 
numerical factors due to a different choice of trial func- 
tions. 

Increasing the Ra number beyond Ra°'^ the oscillat- 
ing frequency lo starts to decrease, until it vanishes at a 
certain Ra2, with the result that for Ra > Ra2 the lin- 
ear instability is stationary again. Ra2 is a complicated 
function of B with Ra2{Bi) = Raf w Raf{Bi). This 
disappearance of the oscillation frequency does not have 
an analog in the case of negative--;/' colloids (without sed- 
imentation) and again, the reason is the independence of 
the sedimentation current from the temperature gradi- 
ent. 

The fully developed state with a concentration gradient 
due to both, sedimentation and thermodiffusion is diffi- 
cult to realize in experiments. Not only one has to wait 
initially for the sedimentation gradient to develop, also 
after each temperature gradient step one has to wait for 
the thermo-diffusive response to be finished. Neverthe- 
less, this case plays an important role in the interpreta- 
tion of the nonlinear bifurcation scenario, below. The lin- 
ear stability analysis of this case is very similar to that of 



the stratified case discussed before, except for the concen- 
tration gradient, which now reads (in dimensionless form) 
d,C = 1- B/iijRa), rather than d^C = -B/{ipRa) as 
before. This difference can be accounted for by substi- 
tuting B with B - %l:Ra in the Eqs.(ITTl)-([14]) leading to a 
linear instability behavior qualitatively the same as be- 
fore. The threshold of the stationary instability is now 



B 



B 



(15) 



which is reduced by a factor oil/ L (at the same B) due to 
the destabilizing effect of thermodiffusion. On the other 
hand, sedimentation has to be stronger by that factor L, 
in order to be relevant. Thus, the oscillatory instability 
with the threshold 



Ra°f^f = Ra° 



27Pr 



14 + 27Pr(0 + 1 
and the frequency at onset 



■(S-Ba) (16) 



6Pr 



14-^27^^(0 + 1) 



{B - B2) 



(17) 



occurs for B > B2 = i?a°-0. The relevant sedimentation 
strength is now B2 ~ 10"*^. Again, the frequency decreases 
with increasing Ra number and vanishes at a certain 
Ra2j leading to a stationary instability. There is again 
a 'triple' point, i?af ^/(Bs) = Ra^'^^i {B2) = Ra2{B2) and 
for very large B, Ra2{B) (x {B — B2)/{ip + 1), asymptot- 
ically. Eqs. dUD-dni correspond to Eqs. (26) and (27) 
of Shliomis & Smorodin [2005]. 

4. Bifurcation Scenario 

A linear theory can neither predict the actual geometry 
of the emerging pattern, nor can it give the its type and 
its saturation amplitude obtained in the long time limit. 
It also fails to describe the complete bifurcation topol- 
ogy, in particular it misses finite amplitude instabilities 
and sometimes it delivers (linear) instabilities that turn 
out to be transients, only. This happens for the thermal 
instability in colloidal suspensions (without sedimenta- 
tion), in the case of a negative separation ratio ^, where 
a linear oscillatory instability relaxes back to the quies- 
cent stable state [Ryskin & Pleiner 2005] . Since a similar 
linear oscillatory instability has been found in the previ- 
ous section, it is compulsory in the present case to discuss 
the nonlinear behavior. This can be done most easily by 
numerical methods and will be described first, reveal- 
ing in particular the time evolution into the saturation 
state. Afterwards, we present an approximate analyti- 
cal expression of the amplitude of the (stationary) state 
that covers fairly well the important nonlinear features 
and allows for their understanding. 

We start with the same numerical method as already 
used in [Ryskin & Pleiner 2005], section IV-A, to inves- 
tigate the time evolution of the system in the nonlin- 
ear regime. This method is essentially an extension of 
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the Lorenz model [Lorenz 1963, Veronis 1965, Ahlers & 
Liicke 1985]. The simphfying idea of this method is to 
look only for solutions corresponding to a 2-dimensional 
pattern (convection rolls). From the thermal instability 
behavior of colloidal binary mixtures one knows that the 
roll pattern is unstable against a square pattern (and for 
high ip values to certain kinds of cross-rolls) close to Ra^ 
[Huke et al. 2007]. This might be still the case when 
sedimentation is included to the analysis, although no 
studies on that are available. Nevertheless, the basic no- 
tions of the nonlinear instabilities, which we will derive 
below, are certainly also true for square patterns, but 
much harder to get numerically and analytically than for 
the roll pattern. 

For convection rolls, which are periodic with wave 
number k in the lateral direction, the ansatz reads 

C{x,z,t) — Co{z,t) + ci{z,t) coskx, (18) 
T{x,z,t) = ~z + eo{z,t) +0i{z,t)coskx, (19) 
w{x,z,t) = wi{z,t) coskx. (20) 

taking into account five modes. The crucial difference 
to the previous work is that here we have an inhomoge- 
neous boundary condition, Eq. ([7]), for the concentration 
variable. 

The initial state has zero velocity and a linear temper- 
ature profile, while for the concentration field we assume 
consecutively a homogeneous, a stratified, and a fully de- 
veloped field, with dzC — 0, dzC = 1 — B/{'tpRa), and 
dzC = —B/^ipRa), respectively, as discussed in the pre- 
vious Section. We then apply additionally a small per- 
turbation of the velocity field wi of the form ~ cos^(7r2:). 
With typical parameter values, ijj = 10, Pr = 7, and 
L = 10^**, and for either moderate B = 400 or rather 
high values B — 18000 a Ra number is chosen above 
the linear threshold value, thus monitoring the tempo- 
ral development of the instabilities. In almost all cases 
we find that the system approaches a stationary state at 
long times, independent of the initial state (a rare excep- 
tion is discussed below). However, how this stationary 
state is reached depends strongly on the linear behavior: 
If for B = 400 the Ra number is chosen to be in the 
linear oscillatory regime, between i?a°* < Ra < Ra2, 
Eq. p^ . the amphtude first starts to oscillate around 
zero, before it increases and, after some overshoot wig- 
gles, reaches the asymptotic constant value, cf. FiglT^). 
Choosing instead a Ra number above Ra2, where the 
system shows a stationary instability, linearly, the am- 
plitude increases directly from zero to its final values, cf. 
FigHb). This explains the reduction of the oscillation 
frequency with increasing Rayleigh number (Fig. 1 of 
Shliomis & Smorodin [2005] ) as a remnant of the linear 
behavior rather than a genuine nonlinear effect. In addi- 
tion, we show that this transient frequency becomes zero, 
when the Rayleigh number exceeds Ra2- 

Similarly, and even more pronounced is this behav- 
ior for large values B > B2 and a Ra number above 
jjqOs,/ RaP: Starting from the homogeneous ground 



6 


A 






4 
2 


/ 

/ 

/ / 




t 


-2 


2 


4 6 8 


10 


-4 






a) 





/ 

/ / 






8 
6 


1 / 
/ / 






4 


/ / 






2 


/ / 




h) 




0.5 


1.0 1.5 


t 



FIG. 1: The time evolution of the convection amplitude for 
B — 400 at two different values of the Rayleigh number corre- 
sponding a) to the oscillatory linear regime, Ra — 2200, and 
b) to the stationary linear regime, Ra = 3200. The broken 
lines are appropriate reference solutions for the case without 
sedimentation {B = 0). 



state, where the linear stability predicts a stationary in- 
stability, since i/j is positive [Ryskin et al. 2003], the 
nonlinear numerical solution shows a smooth and mono- 
tonic transition to the final value of the amplitude, cf. 
FigHb). Starting on the other hand from the fully devel- 
oped ground state, where linear theory gives an oscilla- 
tory instability, Egs. lfTB)) and p7|) . the (nonlinear) ampli- 
tude oscillates for a rather long time around a zero value 
with growing peak amplitude, until finally it switches to 
the stationary finite value, cf. Fig[2^). Starting from the 
stratified ground state, the system is linearly stable, since 
for such a large B value Ra°'^ ^ i?a°. However, with time 
the concentration profile evolves towards the fully devel- 
oped one giving rise, at the end, to the behavior described 
before. The final state is the same stationary convection 
in all three cases. This behavior is in marked contrast to 
the case of the linear oscillatory instability when heat- 
ing a system with negative separation ratio from below 
(without sedimentation), where the initial convective os- 
cillations relaxes back to the non- convective state [Ryskin 
& Pleiner, 2005]. 

For the convection amplitude of the ultimate asymp- 
totic stationary state we derive an approximate analyti- 
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FIG. 2: The amplitude of a finally steady convection as a func- 
tion of time for different initial concentration profiles: a) the 
fully developed profile, and b) the homogenous concentration 
profile. Parameters are the same in both cases Ra = 1796, 
B = 18000, V = 10, L = 10"*. 



cal expression in terms of Ra and B. Substituting Eqs. 
(fT8)) - (|20p into the nonlinear equations of motion (l2])-(|4]) 
and sorting out the different lateral dependencies for the 
stationary state yields the following system of equations 





= Ra k^idi - 


- ^Ci), 


(21) 


(«7lCl) 


= Ldl(Co + ( 




(22) 






)(ci+0i), 


(23) 


{wi0i) + wi 






(24) 






1, 


(25) 



Equation ((22|) can be integrated once. Taking into ac- 
count the boundary condition Eq. ([7]) the concentration 
profile ci(z) is found to be 



Cl 



2L 

Wi 



5.(Co + 0o)-(l 



B 

4iRa 



(26) 



Far from the boundaries Co and ci are proportional to 
L . This follows from the requirement of Eq. ((23)) to be 
consistent with Eq. ([26|l and by taking into account that 
far from the boundaries the derivatives of the functions 
are small. Thus, in Eq. (pS)) Cq can be neglected except 




FIG. 3: The amplitude of the stationary convection as a func- 
tion of the reduced Rayleigh number e = {Ra — RaP)/RaP 
for different values of the barometric number a) B — > 0, b) 
B = 16000, c) B = B2 = Ro°^ « 17534, A) B ^ 19000. 
Solid lines correspond to stable branches. The direction of 
the hysteresis loop is shown by the arrows. 



close to the boundaries and we get 



Cl 



2L / B 

wi \ ^Ra 



-d 



2 Co 



(27) 



To satisfy the boundary conditions for ci, and to find 
the profile of the concentration field near the boundaries, 
one needs to solve the boundary layer problem. This has 
been done in Ryskin & Pleiner [2004], Appendix A. It 
was shown that the boundary layer depth S ~ L^^^ is 
rather small and its contribution to the amplitude equa- 
tion gives only small corrections L. Therefore, Eq. 
([27| can be used to find the velocity and temperature 
distributions. This is obtained approximately by means 
of the trial functions 

wi = ^ cos^(7rz), 00 = G sin(27rz). 6'i = i^cosTrz. 

(28) 

Substituting those profiles into Eqs. ([21]) , dH]) , and (l25l) 
and projecting these equations onto the weight function 
cos^(7rz), leads to a system of three algebraic equations 
for the amplitudes A, F, G. Solving for A, the saturation 
amplitude of convection, we find the implicit expression 



1 + —A ) {Ra^A - L [i^Ra - B]) = RaA (29) 
5 



relating A = A^ / {?>2n'^) to the driving force Ra, the ma- 
terial parameters i, V"; ^^id B. 

In Fig. [3] the amplitude A is shown as a function of the 
reduced Rayleigh number e = {Ra — Ra^)/Ra^ , where 
Ra^ is the linear threshold for thermal convection of a sin- 
gle component liquid. Two qualitatively different types of 
behavior are found depending on the barometric number 
B as the crucial parameter. First, for B < B2 = Ra'^tp, 
there is a monotonic increase of the amplitude with the 
Rayleigh number to arbitrarily high Ra numbers (lines 



7 




FIG. 4: The oscillatiory nonlinear amplitude as a function of 
time for Ra = 1755 > Ra°, B = 17540 > B2, ip ^ W, and 
L = 10-*. 

a,b,c). These lines start with a vertical slope (not gener- 
ally visible at the scale of the figure) at the threshold £c< , 
i.e. at the linear stability threshold i?a^*'^ = B/tfj (which 
is < Ra'^ in this case). In the opposite case, B > B2 (line 
d), there are two branches of the function A{e) in the in- 
terval £1 < £ < £c>- According to our numerical calcula- 
tions the upper branch appears to be stable and the lower 
one unstable. Again, such lines (in particular their un- 
stable branches) intersect the abscissa vertically at £c>, 
i.e. at the linear threshold Ra^'^'f , which is now > Ra'^. 
This is the scenario of a backward instability, where the 
amplitude takes a finite value Ai at the threshold £1 with 

Al = ^{vm;-i) (30) 

£? « — for l>(3>LiP (31) 
5 

£1 « /3 for /3 < LiP (32) 

where (3 = {B — B2) / B2. Increasing the Rayleigh number 
beyond £1 the actual jump from a convection-free to the 
convective state takes place at a Ra number that is gen- 
erally higher than the Ra number, at which the system 
jumps back from the convective to the convection-free 
state lowering Ra. This hysteretic behavior is indicated 
in Fig. [3] by arrows. 

It may come as a surprise that for vanishing amplitudes 
thresholds are found, which belong to the linear stabil- 
ity analysis of the initially fully developed concentration 
profile (rather than the stratified or homogeneous one). 
The reason is that during the development of the con- 
vective patterns, the slow concentration dynamics also 
evolves until the final state is reached. Thus, the non- 
linear behavior of Fig. [3] is obtained for any initial state, 
even the homogeneous one, but only after a long time. 

The last question to be discussed in this section con- 
cerns the existence of oscillatory, or more general non- 
stationary, nonlinear convective states. Numerically, in 
almost all cases a stationary instability has been found. 
Our approximative analytical solution, Eq. (|29p shows 



that for B < B2 there is always a stationary solution, 
while for B > B2 the finite amplitude stationary state 
only exists for £ ^ £1. Of course, the convection-free 
state A = exists for all numbers Ra and B, although 
it looses its stability at some Rac{B), also depending on 
the initial state, as has been discussed by a linear analy- 
sis above. Now it can happen (only for B > B2) that the 
convection-free state is already unstable with respect to a 
linear oscillatory instability, but no stationary nonlinear 
solution exists, leading to a situation where an oscilla- 
tory instability could be the natural response of the sys- 
tem. This requires that the threshold Ra""'-^, Eq. (|16p . 
for an initially fully developed concentration profile, is 
lower than £1. This is possible in a very narrow param- 
eter range only, where £1 is given by Eq. (|32p . i.e. for 
B — B2 < L1PB2 and Ra — Ra^ < LipRa^. In this param- 
eter range an oscillatory solution has indeed been found 
numerically, which is shown in Fig. [?1 The possibility to 
observe this and other solutions experimentally will be 
discussed in the following section. 

5. Conclusions 

We have shown that effects of sedimentation significantly 
change the linear as well as the non-linear behavior of 
thermal convection in colloidal suspensions. We have 
considered the instability of three possible convection- 
free state - an initially homogenous one, an initially strat- 
ified one and the one with the fully developed concentra- 
tion profile. In the two latter cases the linear instability 
can be oscillatory, if the strength of sedimentation, quan- 
tified by the barometric number, is sufficiently large. The 
nonlinear treatment, however, reveals that the oscilla- 
tions are transient only, finally ending up in a stationary 
convective state, where, for simplicity, we have consid- 
ered roll patterns only. Only in a very narrow window 
in the parameter space non- linear oscillations can exist. 
It is not obvious that it is possible to observe such a 
non-linear oscillatory convection in experiment, since in 
the numerics we had to tune the parameters Ra and B 
with an accuracy up to 0.1% to get this state. On the 
other hand our 5-mode model of the nonlinear evolution 
is certainly not exact and we cannot guarantee that this 
numerical oscillatory solution would also be obtained in 
a more refined model. Finally, even if the oscillatory so- 
lution exists in a real physical system, it is certainly very 
difficult to tune the barometric number in experiments 
with an accuracy of 0.1%, since too many physical effects 
contribute to it. 

In a certain range of the parameter space, in partic- 
ular for very high barometric numbers, the stationary 
convective solution comes in two branches, a stable and 
an unstable one. The former can coexist with the sta- 
ble convection-free state, a situation that leads to a hys- 
teretic behavior in experiments. A hysteretic behavior 
was indeed observed in experiments [Bozhko & Putin 
2003]. Although a direct comparison with our theory is 
not possible, since the barometric number in the exper- 
iment is not really known and most of the experimental 
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results are obtained in the presence of a magnetic field, a 
situation which we have not considered so far. We hope 
that our investigations will motivate further experiments 
on sedimentation effects in the convective instabilities of 
colloidal suspensions. 
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